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ABSTRACT 

Emission from the relativistic jet located at hundreds of kpc from the core of the superluminal quasar PKS 
0637-752 was detected at 3.6 and 5.8 /im with the Infrared Array Camera (IRAC) on the Spitzer Space Tele- 
scope. The unprecedented sensitivity and arcsecond resolution of IRAC allows us to explore the mid-infrared 
emission from kiloparsec-scale quasar jets for the first time. The mid-infrared flux from the jet knots, when 
combined with radio and optical fluxes, confirms a synchrotron origin of the radio-to-optical emission and 
constrains very well the high energy end of the nonthermal electron distribution. Assuming the X-rays are 
produced in the relativistically moving knots via inverse Compton scattering of cosmic microwave background 
(CMB) radiation, the infrared observation puts constraints on the matter content of the quasar extended jet. 
Specifically, pure e^e~ pair jet models are unfavorable based on the lack of an infrared bump associated with 
"bulk Comptonization" of CMB photons by an ultrarelativistic jet. 

Subject headings: galaxies: jets — infrared: galaxies — quasars: individual(PKS 0637-752) — radiation 
mechanisms: non-thermal 
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1. INTRODUCTION 

Since th e discovery o f an X-ray jet of the quasar PKS 
0637-752 JChartas et aljr2 000: Schwartz et al. 2000), a num- 
ber of X-ray jets extending to hundreds of kpc distances from 
the quasar nucleus have been unveiled by Chandra (see e.g., 
ISambruna et al.^2004; Marshafl et al 2005). The X-ray emis- 
sion mechanism, however, still remains unsettled for most of 
the large-scale jets of powerful quasars. Based on the spec- 
tral energy distribution (SED) from radio, optical, and X-ray 
bands, it has been argued that the X-ray intensity is too high 
to be explained by synchrotron or synchrotron self-Compton 
radiation from a single population of electrons. 

A currently favored hypothesis for the strong X-ray 
emission is relativistically-amplified inverse Compton (IC) 
scattering of cosmic microwave background (CMB), in 
which the bulk flow of the jet is assumed to be rel- 
ativistic, with a Lorentz factor F ^ 10 all the way to 
nearly Mpc distances, and to be directed toward the ob- 
server at a small viewing angle of 6 ^ F"' ([Tavecchio et aP 
llOOO; Celotti. GhiselHni. & Chiaberge 1200 ll However, 
some difficulties in this model have been recognized, 
and alternative scenarios have been proposed (Aharonian 
|2p02; Dermer & Atovan 2002; Stawarz & Ostrowski 2002; 
lAtovan & Dermer..2004) . Discussions of these models are 
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hampered by poorly known physical conditions in quasar jets, 
as well as by limited observational windows available so far 
In this Letter, we present Spitzer IRAC imaging at wave- 
lengths 3.6 and 5.8 /im of the jet of the superluminal quasar 
PKS 0637-752. With the unprecedented sensitivity afforded 
by IRAC, the observation aimed at measuring the mid- 
infrared part of the broadband nonthermal spectrum of the jet 
knots. Also we have searched for a possible infrared bump, 
which is expected if the jet is indeed highly relativistic at 
hundreds kpc distances from the central engine, and contains 
cold (a Lorentz factor 7^1) electrons. The redshift of PKS 
0637-752 is z = 0.651 (Savage, Browne, & Bolton 1976), so 
we adopt a luminosity distance of Dl = 1.26 x 10^*^ cm, for 
a ACDM cosmology with = 0.27, JIa = 0.73, and Hq = 
71 kms"' Mpc"'. 

2. OBSERVATIONS AND RESULTS 

We have observed PKS 0637-752 with Spitzer IRAC on 
2005 March 27 as part of our Cycle- 1 General Observer pro- 
gram (Spitzer Program ID 3586). IRAC is equipped with 
a four-channel camera, InSb arrays at 3.6 and 4.5 /im and 
Si: As ar rays at 5.8 a nd 8.0 /im, each with a 5f2 x 5.'2 field 
of view jFazio et alj [2004). Only the pair of 3.6 and 5.8 fim 
arrays, observing the same sky simultaneously, was chosen 
for the observation of PKS 0637-752, to obtain longer expo- 
sures in one pair of bandpasses as opposed to two pairs with 
truncated exposure time. The 3.6/5.8 pair was chosen for a 
better spatial resolution and sensitivity. The pixel size in both 
arrays is ~ 1"22. The point spread functions are l."66 and 
1"88 (FWHM) fo r flie 3.6 and 5.8 /im bands, respectively 
JFazio et all 2004Y A total of 50 frames per IRAC band, each 
with a 30-s frame time, were obtained. The 50 images of the 
Basic Calibrated Data processed in the Spitzer Science Center 
(SSC) were combined into a projected image using the SSC 
software mopex. 

Our aims are to explore possible infrared emission of the jet 
close (< 12") to the quasar core. For this purpose, the wings 
of the bright core must be subtracted carefully. We made use 
of the point response function (PRF) image of a bright star 
provided by the SSC to cancel out the PRF wings. We also 
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Fig. 1. — IRAC images (false color) of the PKS 0637-752 jet at 3.6 /im (left) and 5.8 /jm (right) in units of MJy sr '. An angular scale 1" corresponds to 
6.9 kpc projected size. Superimposed on the IRAC images are ATCA 8.6 GHz radio contours (using a circular restoring beam with Of 88 FWHM), starting at 
0.003 mJy beam"' with an increment of 0.003 mJy beam"' ILovell et all2000t) . The emission east of WK5.7 in the 5.8 /im image is confused by residuals from 
wing cancellation. 



employed three bright field stars in our 3.6 /im image as PRF 
templates, and found that the results from the SSC template 
are consistent with those from the field star templates. In our 
5.8 /im image, on the other hand, there are no suitable stars for 
such comparisons. In what follows, we present both the 3.6 
and 5.8 /im results making use of the SSC template alone. The 
infrared fluxes from the core of PKS 0637-752 were mea- 
sured as 5.0 and 9.6 mJy in the 3.6 and 5.8 /im bands, respec- 
tively. These values are below the saturation limits for a 30-s 
frame. 

Figure [U shows the IRAC images of PKS 0637-752 after 
subtraction of the PRF wings of the core. The pixel size of 
the projected images of both channels was set to be 1/4 of 
1 "22. Regions close to the quasar core are blanked out due to 
residuals in the wing-subtraction. Mid-infrared counterparts 
of the main optical knots, WK7.8 and WK8.9 ( S chwartz et al] 
llpOO), which are also the main X-ray features, are clearly vis- 
ible in both channels, located ^ 8" west of the core. The 
two optical knots are marginally separated in the infrared im- 
ages. In the 3.6 /im image, WK8.9 appears to be brighter than 
WK7.8, as in the optical. In the X-ray, WK8 .9 is brighter 
than WK7.8 but the fluxes are more similar (IChartas et alJ 
|2003). No significant infrared emission can be seen at WK5.7, 
WK9.7 3 times fainter than WK8.9 in the optical), or far- 
ther components after a bending point at ~ 10" west of the 
quasar 

Photometry was performed with 5" diameter apertures en- 
closing both WK7.8 and WK8.9, which finds /a.e = 6.5 ± 
0.5 /iJy and g = 11 ±2 /iJy for 3.6 and 5.8 /im, respec- 
tively. The 1(7 errors include the uncertainties associated 
with PRF removal, adopting 10% of the quasar's wing in- 
tensity at the location being considered. Note that we mea- 
sure the sum of infrared flux from the two optical knots. The 
infrared-to-optical slope** is aio — 1.6, significantly steeper 
than the spectral index at either radio or X-ray wavelengths, 

* We define ajo = -log(/ir//„p)/log(i/ir/i^op), with y-„ = 0.83 X lO'* Hz 
(for 3.6 ; ^m) and i/pp = 4.3 X lO'"* Hz. We take f, = 6.5 ^liy and /op = 
0.48 /tJy ISchwartz et all200(D . 



i.e., = 0.81 ± 0.01 and ax = 0.85 ± 0.08 JChartas et alJ 

Although the separation of the two knots is only ~ 1", we 
can derive a rough flux ratio between WK7.8 and WK8.9 at 
the 3.6 /im band, thanks to the high significance and better 
PRF in this band. The flux from WK8.9 is higher by a factor 
of 1 .5 ± 0.6 than WK7.8, in agreement with a factor of ~ 1 .4 
in the optical. 

We can place upper limits on the IR fluxes from a "quiet" 
region preceding the nonthermal condensations WK7.8-8.9. 
For a white box in Fig.[nthat encloses WK5.7, /3 e <2 /iJy and 
/5.8 <5 /iJy (at a 2<t level) are obtained. The flux uncertainty 
is computed as l\fNa, where is the "noise pixels" taken 
from Fazio et al. (2004) and cr presents the background noise 
in which PRF-removal uncertain ty is also formally included. 
We use these flux limits later in ^3.21 

3. DISCUSSION 

3.1. SEDofKnot WK7.8 

The mid-infrared fluxes that we have measured fill cen- 
tral points in the SED of the jet knots and, therefore, set 
an important constraint on models of the broadband emis- 
sion. In Fig. 12] we show the SED for knot WK7.8, the first 
bright knot which has been frequently modeled in the litera- 
ture. The knot is unresolved both in the radio image (ATCA) 
and the HST image ( Sch wartz et aUBOOOt) . Also the corre- 
sponding knot is visible in the deconvolved Chandra X-ray 
image (Chartas et al. 2000). The emissions in these bands 
most likely emerge from the same physical volume. (In con- 
trast, knot WK8.9 appears elongated along the jet in the op- 
tical, and so may involve multiple knots.) Setting a flux ratio 
in the infrared between WK7.8 and WK8.9 to be 1 .4 based on 
the optical (see above), we associate flux densities of 2.7 /iJy 
(at 3.6 /im) and 4.6 /iJy (at 5.8 /im) with knot WK7.8. 

The radio-to-optical spectrum is now arguably attributable 
to synchrotron radiation, because extrapolation of the radio 
spectrum can smoothly connect with the infrared and opti- 
cal fluxes assuming normal steepening due to radiative cool- 
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Fig. 2.— Broadband SED for knot WK7.8. IRAC 3.6 and 5.8 ^im fluxes, 
2.7 and 4.6 /iJy, respectively, are shown as open circles. The radio flux at 
8.6 GHz, 44 mJy, is extracted from the ATCA map centered on WK7.8 with 
a radius 0"8, while the 4.8 GHz flux is set to satisfy a, = 0.8. The HST 
optical (0.2 fi]y at 697 nm) and Chandra X-ray data (6.2 nJy at 1 keV) are 
taken from the list of Kataoka & Stawarz ( 20Q5J . The solid and dashed hnes 
represent the synchrotron+IC/CMI3 modeling (see text). 



ing. Specifically, we first assume that the energy distribution 
of the electrons (in the jet co-moving frame) follows A^(7) = 
/r7"*exp(-7/7max) for 7 > 7inin, where K is the normaliza- 
tion, 7 denotes the Lorentz factor of relativistic electrons, and 
s = 2.6 to match the radio slope a ~ 0.8 (s = 2a+ 1). To repro- 
duce the optical flux, we obtain 7inax — 3.6 x 10^((5ioB-5)~'''^, 
where 610 = S/IQ and B-s = B/(10~^G) are the characteristic 
Doppler factor' and co-moving magnetic field strength (the 
solid line in Fig.|3- The low energy cutoff 7n,in is arbitrary as 
long as 7,„,„ <2 x lO^SioBsT^/l 

This simple model for the synchrotron component overpre- 
dicts the mid-infrared flux by a factor of ^ 2. To account 
for this discrepancy, we modified the form of the electron 
distribution to A^(7) = ^^7~'(1 + 7/7br)~' exp(-7/7max) with 
i = 2.6, which is a smooth broken power law with an ex- 
ponential cutoff, taking account of possible synchrotron and 
IC (in the Thomson regime) cooling. Fitting the SED with 
this model, we obtain 7br ~ 1.4 x 10^ {S\qB-s)~^^^ and 7n]ax — 
1.1 X 10''((5ioB-5)"'''^ (the dashed line in Fig.Q. 

Here we also model the X-ray spectrum, for illustr ative pur- 
pose, by the IC/CMB model (Tavecchio et al] l2000h . namely 
by the model invoking relativistically-amplified IC scattering 
on CMB radiation, with th e Do ppler factor of (5 = 8.4 [cor- 
responding to case (b) in ^3.21 and jj^iin = 20. In order to 
produce the X-ray flux and at the same time to suppress the 
optical flux by the IC process, we obtain a tight constraint, 
10 < 7inin < 40. Note that the IC/CMB model for the X- 
rays generally requires large power carried by nonthermal rel- 
ativistic electrons, Lk^ot — Trr^T^/ScUe, where r = 1 kpc is the 
adopted radius of the knot and is the energy density of non- 
thermal electrons in the jet frame, which can be written as 
Ue = {-f)meC-ne = [(/?- 1)/(/?-2)] 7nim?McC- with the number 
density of nonthermal electrons . For example, one needs 
Lkn ot ^ 9.2 X lO'*^ ergs s"' in the case of the parameter set (b) 
in ^^{d = 8.4, r = 12, and 7n,in = 20). 

The detected mid-infrared fluxes well constrain the highest 
energy part of the electron distribution, roughly 7 10^~^. In 
addition, as we show below, the non-detection of infrared light 
from the jet preceding WK7.8 strongly constrains the lowest 

' The Doppler factor is defined as <5 = [r(l — /3cosS)]"' where jic the 
velocity of the jet, V = (\—0^'y^l^ the bulk Lorentz factor of the jet, and is 
the observing angle with respect to the jet direction. 



energy (7 ~ 1) electrons, which have yet to be accelerated. 

3.2. Bulk Compton Radiation 

The relativistically-amplified IC/CMB scenario raises the 
interesting possibility that an infrared bump may appear as 
a result of "bulk Comptonization" (hereafter BC), namely 
Comptonization of the CMB photons by "cold" (7 1) 
electrons o r positrons in the ultrarelativistic jet stream 
(IsikoraeTal. 1997; Georganopoulos et al. 2005). The BC 
bump is expected to have a blackbody-like spectrum peak- 
ing at the energy £bc — 2(5recMB, independent of redshift, 
where Ecmb is the mean energy of the CMB radiation at z = 0. 
For r = (5 = 10, the bump appears at the infrared wavelength 
Abc ^10 /im. The isotropic BC luminosity from relativistic 
jet containing cold electrons, assuming a cylindrical geometry 
of radius r and length ZcoW, is given bv lGeor^anopoulos et alJ 
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where aj is the Thomson cross section, Ucmb = 4.2 x 
10~'^(1 +z)* ergs cm"-', and Lcoid is the jet power carried by 
cold electrons that upscatter CMB radiation into the infrared, 
defined as Lcoid = 7rr^r^/3c«coidWJeC^ with the number density 
of cold electrons in the jet frame, Wcoid- 

Along a "quiet" stream between the core and the first bright 
knot WK7.8, it is expected that unseen cold particles are prop- 
agating and generating BC radiation. In Fig. |3l we present 
BC spectra calculated for the quiet portion of the jet, assum- 
ing for illustrative purpose that Lcoid = lO'*^ ergs s"'. The 
length of this region is set to be /coid = (13.8/ sin 6*) kpc ap- 
propriate for the white box in Fig. ^ where the flux upper 
limits were derived in ^ We consider the three charac- 
teristic cases as foflows: (a) T = 12, 9 = 3.5°, (b) T = 12, 
9 = 6.5°, (c) r = 15, 6* = 8.3°. These parameters ai-e cho- 
sen for the consistency with the VSOP and VLBI observa- 
tions, which found superluminal motions on milliarcsecond 
scale (~ 100 pc distances from the core) with a m ean apparent 
speed of /3appC = 1 1 .4±0.6c, indicating T > 1 1 .4 jLovell et all 
mmh . The viewing angles are set by the well-known re- 
lation (seelUrrv & Padovani 1995) /3app = ^9 sin 61/(1-/3 cos 6*) 
with /3app = 11.4, which yields two solutions for given F. For 
r = 15, the smaller angle solution 9=l.E° requires an unlikely 
long jet, exceeding a de-projected total length of 2.2 Mpc, and 
therefore is omitted. From Fig.|3] one can constrain Lcoid by 
comparing the calculated BC spectra to the flux upper limits 
obtained from the IRAC data. The results are summarized in 
TablelU 

In Tabled the jet power in nonthermal electrons, L^not, is 
given based on the IC/CMB modeling of knot WK7.8 (with 
7min = 20), together with the lower limits on ^ = Lknot/^-coid = 
{'i)ne/n^o\d- The ratio of the energy densities of nonthermal 
electrons and magnetic fields ub is also shown. It appears 
that case (c) has unc omfortably large Lknot an d Ue/us because 
of small 5 (see also iDermer & Atovanl2004l) . 

The limits on Lcoid impose meaningful constraints on the jet 
models. If the jet is dynamically dominated by e^e~ pairs, one 
can assume that the energy flux of nonthermal electrons at the 
knot is stored in cold e^e~ pairs before the first knot, therefore 
1^ < 1 (i.e. Lcoid > iknot)- However, we obtain ^ > 1 for all 
cases (see Table Our estimates of ^ cannot be noticeably 
reduced by increasing 7,ni„ (through the relation L^not « 7"°,,^ ) 
as 7min ^ 40 (5^0 is set by the X-ray spectrum. Also, changing 
61 in a plausible range 3° < 6* < 9° both for T = 12 and 15 does 
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Fig. 3. — Bulk Comptonizadon spectra expected from the WK5.7 region 
(2" long) assuming L^oid = lO'** ergs s"' and (a) T = 12, = 3.5° (solid 
curve), (b) T = 12, 6* = 6.5° (dashed curve), (c) T = 15, 6 = 8.3° (dotted 
curve). Also shown is a synchi'otron spectrum (solid horizontal line) extrap- 
olated from the radio flux at 8.6 GHz with the spectral form of the dotted line 
in Fig.m The upper limits on IRAC 3.6 and 5.8 {im fluxes are shown as open 
circles. The op tical upper limit is obtained from the 696.9 nm HST image 
ISchwartz et"ali,20Qffi . 



TABLE 1 

Jet powers in cold and nonthermal electrons. 



Case r (5 Lcoid iimot C "cow/nf «f/«s 

(deg) (ergss"') (ergss"') 

(a) 12 3.5 15.6 < 1.0 X 10''5 2.8 X 10"^ > 2.8 <20 (169 

(b) 12 6.5 8.4 <1.6x 10''M.2x lO*'' >5.8 < 10 92 

(c) 15 8.3 5.3 < 1.4x10''^ 1.9x10"'' > 14 <4 2.9 x 10^ 



that jet composition with pure e^e~ pairs can be excluded by 
the absence of BC radiation at X-ray energies due to Comp- 
tonization of UV photons. The present work is the first time 
that jet composition on scales of many kiloparsecs has been 
constrained in this way. 

The number density of cold electrons relative to that of 
nonthermal electrons can be constrained as ncoid/we < 4-20 
(see Tabled, through the relation ^ = (7)n(./«coid and (7) ~ 
2.77tnin ^ 50. We note that these limits are difficult to explain 
within the standard picture of particle acceleration, namely 
diffusive shock acceleration, in which it is generally expected 
that only a small fraction of the cold electrons can be accel- 
erated to nonthermal energies, and consequently ncoid/«f> ^ 1 ■ 
One may need to invoke alternative models regarding the X- 
ray emission mechanism (namely IC/CMB) and/or the accel- 
eration mechanism. 

A possible caveat is that the above discussion relies on the 
simplified picture that the power transported by the jet is con- 
stant. If the ejection of power at the base of the jet is discon- 
tinuous, the emission knots may represent "power peaks" and 
the average jet power could be smaller than the value of L^not 
we calculated. It should also be noted that in our calculations 
of the BC spectra, we have assumed monoenergetic cold elec- 
trons of 7 = 1 . If instead a certain fraction of pairs is "hot" 
with 7 of a few, the BC spectrum can be different from our 
calculation, extending to optical energies. Even so, because 
of the lower limit on the optical flux, the limits on Lcoid would 
be similar to, or tighter than, what we obtained above. 



not affect our conclusion in this regard. Therefore, as long as 
the IC/CMB scenario is correct, our result argues against pure 
e^e~ models. 

Interestingly, on sub-pc scales in radio loud quasars, a simi- 
lar conclusion has been reached bv lSikora & Madeiskilll2()()( ?) 
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